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ABSTRACT 



Mid-infrared images frequently suffer artifacts and extended point spread 
functions (PSFs). We investigate the characteristics of the artifacts and the 
PSFs in images obtained with the Infrared Camera (IRC) onboard AKARI at 
four mid-infrared bands of the S7 (7/im), Sll (11 /im), L15 (15/im), and L24 
(24 /im). Removal of the artifacts significantly improves the reliability of the ref- 
erence data for flat-fielding at the L15 and L24 bands. A set of models of the IRC 
PSFs is also constructed from on-orbit data. These PSFs have extended com- 
ponents that come from diffraction and scattering within the detector arrays. 
We estimate the aperture correction factors for point sources and the surface 
brightness correction factors for diffuse sources. We conclude that the surface 
brightness correction factors range from 0.95 to 0.8, taking account of the ex- 
tended component of the PSFs. To correct for the extended PSF effects for the 
study of faint structures, we also develop an image reconstruction method, which 
consists of the deconvolution with the PSF and the convolution with an appro- 
priate Gaussian. The appropriate removal of the artifacts, improved fiat-fielding, 
and image reconstruction with the extended PSFs enable us to investigate de- 
tailed structures of extended sources in IRC mid-infrared images. 

Subject headings: Data Analysis and Techniques 
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Introduction 



Mid-infrared mu 



(e.g. 


Meixner 


2008; 


Ita et al. 


2008) 


2008; 


Onaka et al. 


20101). However, i 



tiband imagin g data offer a useful tool to study dusty stellar objects 



20081 1 as well as the interstellar medium (e.g., 



Bernard et al. 



20101 ). However, mid-infrared observations frequentl y suffer from variou s 



ler et a. 



2004) 



20061), 



instrumental effects, such as the anomalies in the detector performance ( iPip. 
the extended point spread functions (PSFs) due to the diffraction of light (jjarrett 
and artifacts. The corrections for these annoying effects should become crucial particularly 
when we are interested in the faint diffuse objects adjacent to bright objects. In this Paper 
we report detailed investigations on the qualities of the mid-infrared imaging data taken 
with the Infrared Camera on board Japanese AKARI satellite to evaluate their artifacts 
and to offer useful information needed for the advanced data reduction procedures to users 
of the AKARI mid-infrared data, which are now publically available from the data archive 
server. 



The infrared satellite AKARI carried out about 5000 pointed observations in addition 



to the all sky survey at t 



le mi d- to far-infrared bands durin g its cold phase (u ntil 2008 



Aug 24) f lMurakami et al. 



Onaka et al. 



20071 ) on board 



20071 ). The Infrared Camera (IRC: 
AKARI made near- to mid-infrared (2-26 ^m) imaging as well as spectroscopic observations 
in the pointed observation mode. The IRC has three channels: NIR, MIR-S and MIR-L, 
which cover the spectral ranges 1.8-5.5 /im, 4.6-13.4 /xm, and 12.6-26.5 /xm, respectively. 
Each channel has a field-of-view of about 10^ x 10^ and is equipped with three medium-band 



filters for imaging observations (see 



Onaka et al. 



20071 for details). 



When a b right object comes i nto the field-of-view, several artifacts, often termed as 



"ghosts" (e.g. 



Lorente et al. 



20081 ). appear in MIR images of S7 (7/im), Sll (11/im), 



15 (15)um), and L 24 (24 fim) bands due to reflections among the optical elements 



Lorente et al 



20081 ). The scattering of light internal to the Si: As detector arrays of the 
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AKARI also spreads out the signals from the object over the entire array. This scattering 
together with the diffraction produces very extended components in the PSFs. It is not 
straightforward to separate artifacts from the PSF unambiguously. In this Paper, we 
define artifacts as components that move differently from the real object in the image or 
those showing very distinct, asymmetric patterns (§ 12. ip . In our definition, anomalies of 
the detector array and the effects of diffraction are included in the PSF. We present high 
dynamic range (HDR) PSFs, aperture correction factors, and surface brightness correction 
factors for diffuse sources at S7, Sll, L15, and L24 bands (§ 12. 2p . The obtained HDR PSFs 
are useful in the process of image reconstruction (see §12]), which enable us not just to carry 
out an accurate photometry of point sources but also to get much improved information 
on the structures and surface brightness of diffuse extended sources. We also obtain the 
revised reference data for flat-fielding (hereafter fiat data) for L15 and L24 that free from 
the effects of the artifacts (see § 12. 3p . 

Among the information obtained in this study, the information commonly used in the 
process of data reduction (e.g., new flat dataset for L15 and L24) is included in the latest 
update of the imaging toolkit ver. 20110225. On the other hands, further advanced data 
reduction procedures, beyond the treatments made with the standard pipeline, may be 
necessary for the data of relatively faint extended sources and for those affected by artifacts 
due to an inclusion of bright objects in the FOV. This paper aims to provide the information 
needed in the course of such advanced data reduction procedures to the general users of the 
AKARI MIR datasets and to demonstrate the method of imaging reconstruction and its 
application to the AKARI MIR data of M81 to exemplify the usage of information newly 
given in this paper. 
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2. Detailed examination on the image qualities of AKARI/IHC data 



2.1. Artifacts in IRC images 



2.1.1. MIR-S artifacts 



In the MIR-S channel, a point-like artifact appears with an interval of about 24 pixels 
from the corresponding real point source in the scan direction (hereafter, defined as the 
Y-axis), which is attributed to internal refiections in the beam splitter since it does not 
split the light perfectly (Figured]). We employ archival data of CRL618 and Red Rectangle 
Nebula of Open Time Programs as a calibration source for the artifact, which are listed in 
Table [H The outskirts of the PSF are fairly symmetric about the central source and the 
artifact shows a distinct, point-like pattern. Using this characteristic, we derive the artifact 
patterns at the S7 and Sll bands by subtracting the signals 24 pixels away from the real 
source in the opposite side from the image as the background. No significant difference is 
noticed in the spatial pattern of the artifact between the two objects. 

Despite the consistency in the spatial pattern of the artifact, the ratio of the intensity 
of the artifact to that of the original object varies among objects. Since the degree of 
imperfection va ries with wavelen gth, the intensity of the artifacts can have strong spectral 



dependence (cf. 



Okumura 



20031 ). We investigate the wavelength dependence of the artifact 



by using spectroscopic 



emission ( lOnaka et al. 



datasets. The IRC has a n arrow slit (Ns) for spectroscopy of diffuse 



2007 



Ohyama et al. 



20071 ). The dispersion direction is along the 



Y-axis and the artifact of an emission line appears at a different wavelength. Three datasets 
of spectroscopic observations of two Galactic H II regions and a planetary nebula are 
employed to investigate the spectral dependence (see Tabled]). We use the lines of [Ar III] 
8.99 /xm, [S IV] 10.51 /im, [Ne II] 12.81 //m, and the unidentified infrared (UIR) band at 
11.3 fim. The relative intensity of the artifact varies with the wavelength of the incident 
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photon from 1 % to 4.5 % for 9 to 13 /im (see Figure [2]). The tendency of the wavelength 
dependence is roughly consiste nt with the wavele ngth dependence of transmission of 
the beam splitter (Figure 5 of 



Onakaetal 



20071 ). Because of this characteristic, the 



intensity of the artifact depends on the spectrum of the object and cannot be estimated 
straightforwardly. In practice, the estimate of the intensity of this artifact has to be made 
interactively. Examples of the relative intensity of the artifacts for various objects are given 
in Table |2l The contributions of the artifacts to the total input signals are measured to 
be less than 1 % at S7 and about 1-3 % at Sll, the latter of which is consistent with the 
values obtained from the spectroscopic data. 



2.1.2. MIR-L artifacts 

MIR-L images are severely affected by artifacts originating in reflections among the 
optical elements and the surface of the detector array because of the imperfect AR coating. 
There are at least three types of artifacts (A,B and C) overlapped with the extended PSF 
(§4) in MIR-L images (Figure |3]), one relatively compact (artifact B) and the others quite 
extended (artifacts A and C). The optical path of each artifact is identified by ray-tracing 
calculations. Their major paths are the reflection between the surface of the detector array 
and the KRS-5 lenses, whose AR coating is not perfect particularly at the longer wavelength 
edge of the spectral range of the MIR-L. Thus the artifacts appear more strongly at L24 
than at L15. The artifact images from a point source are not focused onto a point, but 
appear as a diffuse pattern. We examine the properties of these artifacts and derive the 
surface brightnesses and the positions relative to the original object using the on-orbit data 
listed in Table [T] Artifacts A and C, and the PSFs (see § 12.2. ip are very extended and 
the subtraction of the background sky needs to be made with high accuracy. There is no 
appropriate position for the estimate of the background sky in the image, and we use the 
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data taken at a near sky position that do not have bright objects in their field-of-view as 
the sky background. 



20071), a 



During an pointed observation carried out with AOT IRC02 (lOnaka et al. 
dithering operation was carried out to settle the target at three different positions on the 
detector array. Using these dithered data, the movement of the artifacts is estimated by a 
least-square fit for the region where each component of the artifacts does not overlap with 
each other. The relative movement of each artifact is assumed to be proportional to the 
movement of the original object. Thus the relative movement of the reference position of 
each artifact {xr,yr) from the origin {xo,yo) is given by the position of the incident object 

{Xr - Xo, Vr - Uo) = {a{xi - Xo),b{yi - Uo)), (1) 

where a and b are the fitting parameters. We set the origin as the center of the array, 
{xoiVo) = (128 pix, 128 pix). The artifact patterns do not change with the position of the 
object in the field-of-view and thus the reference position {xr,yr) is defined as a position 
near the center of the artifact. For real objects and their PSFs in the image, a and b should 
be unity, but these are not for the artifacts. By the difference in the relative movement, the 
artifacts can be separated from the PSF pattern. The derived fitting parameters are listed 
in Table [31 

The artifacts on L15 images are too faint to derive reliable position movements, and 
the fit is made only for L24 images. We assume the same values of a and b for L15 as 
derived for L24 since the positions are determined by internal reflections and the spectral 
dependence is thought to be negligible. The parameters a and b for artifact A are both 
positive and thus artifact A moves in the same direction as the real object. Artifacts B 
and C move, on the other hand, in the opposite direction relative to the real object with 
the same amount on the image. Taking account of the difference in the movement on the 
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image, we determine the patterns of artifacts A, B, and C, and the PSFs both at L15 and 
L24 (Figures m and E]). The spider patterns and the banding structures from the detectors 
are included in the PSFs (see Figure [TT]). Because they move with the same amount in the 
same direction, artifacts B and C cannot be separated uniquely from the relative movement. 
The relative intensities between these two artifacts, however, vary between L15 and L24 
(see Table HI). Each artifact image is shifted and stacked by median average, excluding the 
maximum value to remove the other artifacts from the output image. Artifact A in the L15 
image is very diffuse and fully overlaps with the extended component of the PSF, which 
makes it difficult to extract. We assume the same pattern at L24 as for the low frequency 
(> 25 pixels) of artifact A at L15. Artifact C is much fainter than artifact A at L24 and we 
assume the same pattern at L24 as at L15. 

The characteristics of artifacts A, B, and C are summarized in Table HI The total 
contributions of the artifacts to the input signals are 7.3 % and 16.1 % at L15 and L24, 
respectively. The spectral dependence of the intensity of the artifacts relative to the original 
object cannot be investigated with accuracy because of the lack of sufficient data. We 
assume fixed relative intensities of the artifacts as given in Table HI where the artifact 
signals are normalized by the signal of the target point source within a radius of 7.5 pixels 
(see below). Among the MIR-L images we have investigated, the spectral dependence of 
the relative intensity of the artifacts is less than 15 %. Corrections for the artifacts have 
to be made for each dithered image since the artifacts move differently from real objects. 
Figures [6] and [7] show examples of the original images with the artifacts and those of the 
artifacts removed for the carbon star U Antliae at L24, respectively (pointing ID:1710071.1, 



Arimatsu et al. 



20111 ). They show that the artifacts seen in Figure [6] are subtracted almost 
perfectly after the removal process (Figure [7]), demonstrating the applicability of the artifact 
patterns at L15 and L24 derived in the present study. 
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2.2. Extended PSFs of IRC images 



2.2.1. Extended PSFs 



In general, the PSFs of mid-infrared instruments te nd to be extende d due to the 



2004 ). We determine 



diffraction and the scattering within the detector arrays (iHola et aLl 
the PSFs over the entire field-of-view of the four imaging bands (S7, Sll, L15, L24) by 
subtracting the artifacts from cahbration and science observations of stars. The observations 
used to determine the PSFs are hsted in Table [1] To obtain HDR PSFs, we use targets 
with a wide range of brightness. We also employ short-exposure and long-exposure images, 
the latter of which was taken with 28 times longer exposure time than the former. Images 
after the artifact subtraction at each dithering position are shifted and combined with 
an accuracy of 0.1 pixels to produce a brightness map, and we determine the extended 
component of the PSFs down to 10^^ of the total intensity. Figure [TTl shows the PSFs, 
which indicates significant morphological differences among the imaging bands. Especially 
in the PSFs at S7 and L15, a banding structure along the cross-scan direction (hereafter, 
defined as the X-direction) is seen in addition to the structures produced by the secondary 
mirror spiders separated by 120 degrees. A similar banding pattern is seen in Syitzerll llAC: 



images, which is attributed to light scattering inside the detector ( iPipher et al 



2ooi . 



The banding shows a strong spectral dependence and the morphological differences can be 
attributed partly to the internal reflections in the Si:As detector and/or to diffraction o: ' 



light. Uniform offsets that appear under high incoming photon conditions ( jPipher et al. 



20041 ) are not recognized in the IRC PSF object images, when they are compared with the 



reference sky images (§ I2.1.2p . 

The full width of half maximums (FWHMs) of the obtained PSFs are listed in Table |5l 
The obtained valu es of the FWHMs are smaller than the published values (Table 2 of 



Onaka et al. 



20071 ) at the short wavelength channels. The present values are consistent with 



- 10 - 



20071 ). The published values 



the diffraction-limited performance at 7.3 //m f lKaneda et al.l 
are derived from a number of stellar images and the averaging process could increase the 
FWHM. In another word, they provide "average" FWHMs when the standard reduction 
toolkit is used. On the other hand, the PSFs in the present study are derived from the data 
that are carefully shifted and stacked. This process requires for the data with sufficient 
signal-to-noise ratios. When the present HDR PSF is applied to a standard processed 
image, the inner part of the new PSFs needs to be adjusted to the individual image, which 
might be degraded by a less accurate stacking process, while the outer part provides proper 
information on the extended component of the PSFs. 



2.2.2. Aperture Corrections 

The extended component of the PSFs can be characterized quantitatively by the 
encircled energy. It is also useful for practical applications to calculate the aperture 
correction factor AF, which is defined as the reciprocal of the encircled energy normalized at 
a given radius. We set the reference radius as 7.5 pixels, with whi ch the absolute calibration 



of the MIR-S and MIR-L data is carried out (ITanabe et al. 



20081 ). This radius corresponds 



to 17". 5 5 for t he MIR-S and 17". 88 for the MIR-L channel, respectively. According to 



Jarrett 



(j2006[ ). we approximate AF as a function of the radius r from the center of the 



object as 



AF{r) = Aexp(-r^) + a 



(2) 



where r is in pixels, and A, B and C are the fitting parameters. They are given in Table E] 
and AF{r) is plotted in Figure [T2l Figure [T2] indicates that equation (|2]) provides good fits 
for r > 7.5 pixels. 



The calibration for diffuse emission needs to take account of the fact that a large 
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fraction of the incident photons spread over the detector array. The correction for 
diffuse emission can be given by equation ([2]) with r oo. Thus parameter C gives the 
correction factor for diffuse emission. As indicated in Table O the correction factor has 
spectral dependence. Scattering within the detector array is kn own to have strong spectral 



dependence, peaking in the 5-6 /xm region (jPipher et al. 



20041 ). It explains part of the 



dependence in the MIR-S bands. For the MIR-L bands, diffraction of light is likely to make 
a significant contribution. 



2.3. Revised Flat Frames for MIR-L images 

Diffuse emission images are affected not only by the artifacts and the extended PSF 
components produced by bright objects in the field-of-view, but also by those produced 
by diffuse sky background. Signals of the artifacts produced by the diffuse background in 
each pixel are small, but the total amount of the signals integrated over the entire array 
becomes significant. Flat data for mid- infrared instruments are usually produced from the 
sky background data. Previous flat data for the MIR-L channel are created from a large 
number of background sky images. Thus the artifacts at L15 and L24 can affect the derived 
flat data. In fact, th ey show strange patterns, which are attributed to the artifact effects 
(ILorente et al.ll2008l ). The flat data show a steep decline at the edges of the fleld-of-view 
(Figure IHk) • It is most likely that the decline is produced by the effect of the artifacts, 
especially artifact A. 

To investigate the effects of the artifacts on the flat data, we carried out observations 
of the MIR-L standard star IRAS F06009— 6636 at 25 different positions of a 5x5 grid 
separated by ~ 2' in the fleld-of-view of the MIR-L (Table [1]) and compare the variation in 
the signal with the flat data. The flux of the star is measured with aperture photometry 
with a radius of 7". 5 and is derived for each dithered image. The comparison indicates that 
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the previous flat data overcorrect the flux of a point source at an edge of tlie field-of-view 
by 15 % at most at L24 (Figures IHb andfTOb). Tlie standard deviation of the flux variation 
increases from 3.8 % (without flat-fielding) to 6.4 % if the flat correction is made. These 
observations were useful to estimate the accuracy of the previous flat-fielding process of 
MIR-L images, but did not have sufficient spatial information to produce new flat data. 

We produce new flat data for L15 and L24 from the previous flat data by correcting 
for the artifacts derived in § 12.1.21 First we assume that the sky background is uniform 
over the entire field-of-view and integrate the contribution from the signal of each position 
to obtain the artifact pattern that is produced by the uniform background. Next the 
artifact pattern is searched for in the previous flat data and the amplitude of the artifact 
is estimated. The scaled artifact pattern is then subtracted from the previous flat data. 
The artifact-subtracted-flat data are renormalized at the center position to keep the present 
absolute calibration valid since most of the standard stars are observed around the center 



of the field-of-view (ITanabe et al 



20081 ). Finally we apply the newly-derived flat data 
to the 5x5 standard star observations. The results indicate that the flat is significantly 
improved in the Y-direction (Figure [TOb ) , but there still remains a shallow decline around 
the edge in the X-direction (Figure [9t). There may still be a very diffuse artifact that 
remains uncorrected. We fit the pattern in the X-direction with a quadratic function of the 
position and apply it to the fiat data to create the final fiat data (Figure |8]i). With the 
final artifact-subtracted and slope-corrected flat data, the standard deviation of the flux 
variation of the standard star decreases from 3.8 % to 2.8 % for L24 (see Figures [9li and 

For the L15 band, a similar trend is seen in the previous flat data. The same procedure 
is applied as for the L24 band to obtain new flat data after the slope correction in the 
X-direction. The standard deviation of the flux variation of the standard star improves 
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from 3.2 % (without flat-fielding) to 1.2 % (with the artifacts-subtracted fiat data with the 
slope correction) for the L15 band. 



3. Data Reduction Techniques beyond the Imaging Pipeline Procedures for 

the Advanced Analyses of Diffuse Sources 

A significant fraction of IRC observations in the cold phase were made with the S7, 
Sll, L15, and L24 bands, which provide 
prop erties of the interstellar matter (e.g., 



us useful information to charac 


;erize 


ihe emission 


Onaka et al. 


2009; 


Sakon et al. 


2007; 


Onaka et al. 



20101 ). However, the differences in the PSF sizes as well as in the distinct patterns in a 
large scale, as shown in Table 5 and Figure 11, hamper us from directly making color 
analyses among those bands especially for extended objects such as ISM structures of 
nearby galaxies. Figure [T3] compares the radial profile of the PSF at L24 and that at Sll 
convolved with a Gaussian to have the same FWHM as the L24 PSF (6". 7, see Tabled]). 
It shows that the radial profile of the adjusted PSF of Sll is different from that of L24, 
especially at larger radii beyond 10". In Figure [T4l we exemplify the case of imaging data 
of the nearby galaxy M81. Ring-like artifacts appear around H II regions and the nucleus 
of the galaxy in the L24/S11 flux ratio map of M81 made from this simple convolution of 
the Sll image (Figure [Hb. indicated by the arrows). The ring-like structures come from 
the diffraction pattern at L24, which should be corrected for by making use of the proper 
PSF pattern. To make comparative studies with the four bands it is necessary to adjust the 
PSF pattern properly. We investigate a method to convolve the deconvolved image with an 
appropriate kernel. 



Following 



Gordon et al. 



(120081 ) ■ we define the convolution kernel K{x,y) as 



K(x ,) = DFT-|H!ZP^!:i(M)l 



(3) 
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where DFT indicates the discrete Fourier transform operation and DFT~^ is the reverse 
operation. PSFi is the PSF of the input image and PSF2 is a simple PSF to which all band 
images will be adjusted. The output image Jout(-^, ^) is then obtained by 

Iout{X,Y) = DFT~'\DFT[K{x,y)] ■DFT[h,{x,y)]y (4) 

In these operations, the sampling scale for Ii^{X,Y), /out(-^, ^) and K{x,y) is set as a 
half of the intrinsic pixel scale. The Gaussian profile with the FWHM of 7"is employed for 
PSF2. The radial profile is shown by the blue line in Figure fT3l The convolution with a 
Gaussian suppresses spurious features found in a simple deconvolution operation without 
losing spatial information significantly. 

The present reconstruction procedure provides images that can be used reliably for 
the study of faint extended emission. Figure [14] shows how the reconstruction procedure 
described above improves the image quality taken with the IRC. For example, some discrete 
H II regions are seen as bright point-like spots over the relatively faint diffuse emission 
from the ISM in the M81 in the L24 image in Figure [T^ . Due to the extended component 
of the PSF, especially to the diffraction ring structure of the PSF at L24, the diffuse 
emission around such bright sources are overlapped by the extended artifacts that come 
from the bright sources. Figure 15 shows the surface brightness profiles along a certain 
selected line defined in both images of region A before and after processing the present 
image reconstruction (see Figures 14c and 14d. respectively). A typical level of the diffuse 
emission in region A is about 0.7 MJysr"^ if the area that is not affected by bright sources. 
On the other hand, a bright infrared point-like source centered at Ay = 12.5 with a flux 
density of 0.02 Jy at L24, for example, produces the ring-like pattern that makes a false 
contribution to the surface brightness as bright as 2 MJysr"^ in the profile ranges of 
Ay =7-9 and Ay=16-18 (see Figure 15). In this case, about one third of the whole pixels 
in region A are almost comparably affected by the extended PSF components from the 
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bright sources, which makes it difficult to estimate the emission from the diffuse sources in 
concern. These artifacts are properly removed by the reconstruction procedure, as shown 
in Figures [T4l i and [151 which offers significant improvement in accurate photometry of the 
surface brightness of diffuse sources located even at a closer region to a point source. After 
the present image reconstruction, the surface brightness of the diffuse emission can be 
measured with an accuracy of 0.2 MJysr~^ at L24 if there are no sources brighter than 0.1 
Jy within 5 pixels from the region in question. We note that the present methods do not 
take account of the changes in PSF shapes within the FOV. Thus the accuracy of the present 
image reconstruction may be slightly degraded for diffuse sources around a point source 
located in the corner of the FOV. Figures [T^ and f indicate that the image reconstruction 
significantly improves the quality of the comparative study of four-band images. As shown 
in the reconstructed L24/S11 color map (Figure fT^). the ring-like structures disappear and 
the contrast of the L24/S11 fiux ratio is clearly enhanced compared to the original color 
map (Figure IT^). 

4. Summary 

Mid-infrared images taken with the IRC onboard AKARI suSei artifacts and extended 
PSFs. The artifacts are investigated at four IRC bands, whose contributions to the total 
input signal are measured to be about < 1%, 1 — 3%, 7%, and 16% at S7, Sll, L15, and 
L24, respectively. These artifacts affect the flat-fielding reference data and new flat data 
are derived by removing the effects of the artifacts. The PSFs and the aperture correction 
factors for point and diffuse sources are obtained at S7, Sll, L15, and L24. To adjust 
the spatial resolution at different bands, an image reconstruction method is developed. 
The present investigation enables us to carry out the comparative study of four-band 
mid-infrared images taken with the IRC. The artifacts and extended PSFs of the IRC data 
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have common characteristics inherent in mid-infrared instruments. The present analysis 
will be useful not only for the analysis of IRC MIR data, but also for the design, ground 
tests, on-orbit calibration plans, and data analyses of mid-infrared instruments of future 
facilities. The image reconstruction developed here can also be applied to data taken with 
other instruments. 

The MIR-L ghost patterns, the new flat data of L15 and L24, and the HDR PSFs of 
S7, Sll, L15, and L24 are available at jhttp : / / www. ir . isas . j axa . j p/ AK ARI / Observation / [ 
The new flat data are included in the latest version of the imaging toolkit (ver. 20110225). 
The v alidity of the absolu te calibration is confirmed with the data of standard stars in the 



LMC flTanabe et al. 



20081 ) processed with the latest imaging toolkit with the new flat data. 
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Table 1: List of the datasets used for PSFs, artifacts and flat-fielding investigation 



Observation ID 


AOT" 




Date 


Object 


Note 


MIR-S channel 














4080002.1 


IRC03 


a 


L 


2007 Mar 05 


CRL618 


PSF & SCL 


4080019.1 


IRC03 


a 


L 


2007 Mar 05 


CRL61S 


PSF 


4080007.1 


IRC03 


a 


L 


2006 Sep 28 


Red Rectangle 


PSF & SCL 


4080006.1 


IRC03 


a 


L 


2006 Sep 28 


Red Rectangle 


PSF 


1800623.1 


IRC02 


a 


L 


2007 May 10 


HIP 28909 


PSF 


5020008.1 


IRC03 


a 


L 


2006 Apr 22 


HIP 28909 


PSF 


5124007.1 


IRC03 


a 


L 


2006 Jun 24 


KF01T4 


PSF 


5124019.1 


IRC03 


a 


L 


2006 Aug 03 


KF03T2 


PSF 


1400229.1 


IRC04 b;Ns 


2007 Feb 06 


CAR 057TO005 


SCL 


1401033.1 


IRC04 b;Ns 


2007 Apr 06 


CRU 032 S003 


SCL 


5020048.1 


IRC04 


a;Ns 


2006 Apr 29 


NGC6543 


SCL 


MIR-L channel 














1711365 


1 


IRC02 


a 


L 


2007 May 12 


IRAS 22396-4708 


PSF & SCL 


4080019 


1 


IRC03 


a 


L 


2007 Mar 5 


CRL618 


SCL 


4080002 


1 


IRC03 


a 


L 


2007 Mar 5 


CRL618 


PSF& SCL 


5124105 


1 


IRC03 


a 


L 


2007 Jul 12 


IRAS F06009-6636 


PSF & FLAT 


5124089 


1 


IRC03 


a 


L 


2007 Jun 19 


IRAS F06009-6636 


PSF & FLAT 


5124090 


1 


IRC03 


a 


L 


2007 Jun 22 


IRAS F06009-6636 


PSF & FLAT 


5124091 


1 


IRC03 


a 


L 


2007 Jun 22 


IRAS F06009-6636 


PSF & FLAT 


5124099 


1 


IRC03 


a 


L 


2007 Jul 4 


IRAS F06009-6636 


PSF & FLAT 


5124092 


1 


IRC03 


a 


L 


2007 Jul 21 


IRAS F06009-6636 


PSF & FLAT 


5124098 


1 


IRC03 


a 


L 


2007 Jul 4 


IRAS F06009-6636 


PSF & FLAT 


5124097 


1 


IRC03 


a 


L 


2007 Jul 1 


IRAS F06009-6636 


PSF & FLAT 


5124100 


1 


IRC03 


a 


L 


2007 Jul 5 


IRAS F06009-6636 


PSF & FLAT 


5124083 


1 


IRC03 


a 


L 


2007 Jun 1 


IRAS F06009-6636 


FLAT 


5124081 


1 


IRC03 


a 


L 


2007 Jun 2 


IRAS F06009-6636 


FLAT 


5124082 


1 


IRC03 


a 


L 


2007 Jun 2 


IRAS F06009-6636 


FLAT 


5124084 


1 


IRC03 


a 


L 


2007 Jun 7 


IRAS F06009-6636 


FLAT 


5124085 


1 


IRC03 


a 


L 


2007 Jun 12 


IRAS F06009-6636 


FLAT 


5124087 


1 


IRC03 


a 


L 


2007 Jun 13 


IRAS F06009-6636 


FLAT 


5124086 


1 


IRC03 


a 


L 


2007 Jun 14 


IRAS F06009-6636 


FLAT 


5124088 


1 


IRC03 


a 


L 


2007 Jun 18 


IRAS F06009-6636 


FLAT 


5124094 


1 


IRC03 


a 


L 


2007 Jun 26 


IRAS F06009-6636 


FLAT 


5124095 


1 


IRC()3 


a 


L 


2007 Jun 27 


IRAS F06009-6636 


FLAT 


6124093 


1 


IRC03 


a 


L 


2007 Jun 29 


IRAS F06009-6636 


FLAT 


5124096 


1 


IRC03 


a 


L 


2007 Jun 30 


IRAS F06009-6636 


FLAT 


5124101 


1 


IRC03 


a 


L 


2007 Jul 8 


IRAS F06009-6636 


FLAT 


5124102 


1 


IRC03 


a 


L 


2007 Jul 8 


IRAS F06009-6636 


FLAT 


5124103 


1 


IRC03 


a 


L 


2007 Jul 9 


IRAS F06009-6636 


FLAT 


5124104 


1 


IRC()3 


a 


L 


2007 Jul 12 


IRAS F06009-6636 


FLAT 



"Astronomical Observation Template for the IRC and FIS observations. See ASTRO-F Observer's Manual 



for details of the parameters (http : //www. ir . isas . jaxa. jp/AKARI/Observation/ObsMcin/|. 

^PSF indicates the data that were used to derive the PSF. SCL indicates those used to estimate the scattered 
light contribution (artifacts). FLAT indicates those used to estimate the reference data for flat- fielding at 
L15 and L24. 
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Table 2: Relative intensity of the MIR-S artifact 

Object S7 Sll 

A-typestar^ 0.0065 ± 0.0008 0.016 ±0.002 

K-typestari^ 0.0073 ± 0.0006 0.024 ±0.004 

Zodiacle light^ < 0.0027 0.023 ± 0.008 



°HD 28909 (AO star) 
^KF01T4 (K1.5III star) 
"^Calculated from the slit image. 
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Table 3: Parameters for MIR-L L15 and L24 artifact movements 



Artifact 


a [pixels] 


b [pixels] 


A 


0.364 X (1 ±0.051) 


0.194 X (1 ±0.012) 


B and C 


-1.426 X (1 ±0.0049) 


-1.272 X (1 ±0.0028) 



"The fit parameters in equation ([T|). 



Table 4: Properties of the artifacts at L15 and L24 



Artifact 


Total signal to 


Typical radius 


Surface brightness per pixel 




the object signal 


[pixels] 


relative to the object signal 


L15 








A 


2.5 X 10-2 


100 


1 X 10-6 


B 


1.8 X 10-2 


30 


2 X 10-6 


C 


3.0 X 10-2 


150 


6 X 10-^ 


L24 








A 


6.8 X 10-2 


100 


2.4 X 10-6 


B 


4.9 X 10-2 


30 


1.8 X 10-^ 


C 


4.9 X 10-2 


150 


1 X 10-6 
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Table 5: FWHMs of the PSFs 



IRC band FWHM [arcsec] 
S7 2.9 
Sll 3.3 
L15 4.6 
L24 6.7 



Table 6: Best-fit coefficients of the aperture correction factors 

IRC band A" 
S7 0.776 0.371 0.906 

Sll 0.804 0.428 0.926 

L15 1.671 0.417 0.843 

L24 1.764 0.396 0.817 



"The fit parameters in equation ([2]). 



- 23 - 




Fig. 1.— Example of the MIR-S artifact at Sll for CRL 618 (pointing ID:4080002.1). The 
secondary and tertiary artifacts, which come from the first and secondary artifacts, are also 
seen. The image is displayed in a logarithmic scale. 
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10 11 12 
Wavelength [\im] 



14 



Fig. 2. — Spectral dependence of the artifact in the MIR-S images derived from emission 
hnes in the sht spectroscopy. At 10.5 ([S IV]) and 11.3 (UIR) yum the corresponding artifacts 
are not detected and 3-a upper hmits are shown. 




Fig. 3.— Artifacts at the L24 band. The target object is IRAS 22396-4708. The identified 
artifacts (A, B, and C) are indicated. The image is displayed in a logarithmic scale. The 
entire field-of-view (256x256 pixels) is shown. The black regions in the center are saturated 
pixels. The structures seen in the right edge are the slit mask. The diffraction spikes and 
banding (see § 12.21) can also be seen in the figure and are overlapped on the artifacts, which 
can be separated from the artifacts with the separation procedure (see § 12.1.21) and are finally 
included in the PSF. 
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10-» 5x1 0» 5x1 0" 10- 



10« 5x1 0-» 5x10" 10- 




10' 5x1 0-» 



Fig. 4. — The PSF and the artifacts at L15. The object is located at the center (128 pix, 
128 pix). The arrows in the right two figures indicate the direction of the relative movement. 
All images are displayed in a logarithmic scale with the object signal within 7.5 pixel radius 
normalized as unity. 





5x10-" 10-' 10" 10-' 5x10-= 5x10" 10-' m' 10' 5x1 0-» 5x10" 10-' 10" 10» 5x1 0-» 



Fig. 5. — The PSF and the artifacts at L24. The layout is the same as in Figure |H 
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Fig. 6. — Example of the individual dithered images with the artifacts of U Ant at L24 
(ID:1710071.1). The central region of field-of-view (160x150 pixels) is shown. The color 
scales are in units of ADU per pixel and all figures are displayed in a logarithmic scale. 





















m 








50 100 1 50 200 250 300 350 400 450 



Fig. 7. — Example of the images of U Ant at L24 after the artifacts removal. All figures are 
displayed in a logarithmic scale. 
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Fig. 8. — (a) Previous flat data of L24, (b) artifact A predicted from a uniform background, 
(c) artifact B and C patterns predicted from a uniform sky, and (d) flat data of L24 corrected 
for the artifacts. Tlie entire detector array (256 x 256 pixels) is shown. 
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Fig. 9. — Comparison of tlie aperture photometry of tlie same point source (F06009— 6636) 
at 25 different positions at the L24 band. All the plots are made along the X-axis of the 
image, (a) The photometric data without fiat-fielding, (b) Those corrected with the previous 
flat data, (c) Those corrected with the artifact-subtracted fiat data, (d) Those corrected 
with the artifact-subtracted fiat data with the slope correction shown in Figure c (see text). 
The photometry is all carried out with the standard aperture size (7.5 pixels in radius) and 
the flux is normalized by the average of the fluxes of the data taken at the center position. 
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Fig. 10.— 



Same as Figure |9] except that the plots are made along the Y-axis of the image. 




Fig. 11.— PSFs of AKARI/IRC imaging bands at (a) S7, (b) Sll, (c) L15, and (d) L24. 
All figures are displayed in a logarithmic scale with the total intensity of unity. The contour 
levels are set as 2^""'"^^ x 10^^ of the total intensity. 
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Effective Circular Aperture Radius [arcsec] 



Fig. 12. — Aperture correction factors for the S7 (blue), Sll (green), L15 (orange), and L24 
(red) bands. The hues are the fitting results with equation (|2]). The standard photometry 
radius (7.5 pixels) is shown by the violet and black lines for MIR-S and MIR-L, respectively. 
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Radius [arcsec] 



Fig. 13. — Radial profiles of the PSF at L24 (red line), the PSF at Sll (green line) convolved 
with a Gaussian to adjust the FWHM to 6". 7, and the PSF obtained from the deconvolution- 
convolution process (equation ([3])) for all the four bands (blue line). 
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e. L24/Sl'I (tttfore image reconstruction) 



^ e Mjon A 



b. S11 (after irnage reconstruction) 
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I I d. 1^4 (after image reconstrucfl^) 
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Fig. 14. — Example of the image reconstruction process of the IRC mid-infrared data investigated 
in the present study (M81: pointing ID 5020062.1 and 5020063.1). (a) Artifact-subtracted image of 
Sll, (b) Sll image reconstructed with equation (jj]), (c) artifact-subtracted image of L24, (d) L24 
image reconstructed with equation @, (e) the L24/S11 color map produced from the L24 image (c) 
and the Sll image (a) convolved with a simple Gaussian to adjust the FWHM to 6". 7 (green line 
in Figure [T3]) . and (f) the L24/S11 color map produced from the reconstructed images (b) and (d). 
Arrows indicate bright H II regions and the galactic center of M81, where ring-like artifacts appear 
in (c) and (e). After the image reconstruction procedure, the ring-like artifacts seen in (c) and (e) 
disappear in (d) and (f). The blue line in region A overlaid on (c) and (d) indicates the position of 
the surface brightness profiles shown in Figure [15] and the green boxes indicate the regions used to 
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5 10 15 20 25 

AY [pixels] 

Fig. 15. — Surface brightness profiles at L24 along the line shown in Figures [T^ and d. 
The dashed line represents the profile before image reconstruction (Figure [Ht); the solid 
line shows that after image reconstruction (Figure dUi). The pixel scale is 2". 38, which 
corresponds to about 42 pc assuming the distance of 3.6 Mpc to M81. A bright HII region 
is located at AF = 12.5 pixels and the ring-like pattern appears at Ay ~ 7.5 pixels and 
~ 17.5 pixels. 



